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1 Introduction 


The Ruijsenaars-Schneider-van Diejen (RSvD) models P, |2] are among the most intensively 
studied integrable many particle systems, having numerous relationships with different branches 
of theoretical physics and pure mathematics. They had found applications first in the theory 
of the soliton equations P El m El E], but soon they appeared in the Yang-Mills and the 
Seiberg-Witten theories as well (see e.g. p El El ini E]). Besides these well-known links, 
the RSvD systems and their non-relativistic limits, the Calogero-Moser-Sutherland (CMS) 
systems naEniE], have appeared in the context of random matrix theory, too. Making use of 
the action-angle duality between the different variants of the CMS and the RSvD systems, new 
classes of random matrix ensembles emerged in the literature [T^lfTOlITT] . exhibiting spectacular 
statistical properties. Under the name of classical/quantum duality, it has also been observed 
that the Lax matrices of the CMS and the RSvD models encode the spectra of certain quantum 
spin chains, thereby the purely classical models provide an alternative way to analyze the 
quantum systems, without any reference to the celebrated Bethe Ansatz techniques (for details 
see e.g. [m ESI EHl |2I]). It is also worth mentioning that in the recent papers [22l [23] the 
authors have constructed new integrable tops, closely related to the CMS and the RSvD particle 
systems. Besides the Lax representation of the dynamics, in their studies the associated r- 
matrix structures also turn out to be indispensable. 

The characteristic feature the above exciting new developments all share in common is the 
prominent role played by the Lax matrices of the CMS and the RSvD models. However, all 
these investigations are based on the translational invariant models associated with the An 
root system, exclusively. Apart from the technical difficulties, the probable explanation of this 
state of affair is the very limited knowledge about the Lax representation of the RSvD models 
in association with the non-A„-type root systems. Of course, one can easily construct Lax 
representations for both the C^-type and the RC^-type RSvD models by the Z 2 -folding of the 
A2n-i and the A2n root systems, respectively [21]. However, this trivial approach is only of very 
limited use, since the resulting models contain only a single coupling parameter. Nevertheless, 
working in a symplectic reduction framework, in our papers [25l|26] we succeeded in constructing 
Lax matrices for the rational Cn and the rational BCn RSvD systems with the maximal number 
of independent coupling constants. Motivated by the plethora of potential applications outlined 
above, in this paper we work out the underlying classical r-matrix structures and also provide a 
Lax representation of the dynamics for the rational BCn RSvD model with three independent 
coupling parameters. 

Let us recall that the configuration space of the rational BCn RSvD system is the open 
subset 

c = {A = (Ai,..., A„) e M*" I Ai > ... > A„ > 0} C (1.1) 

that can be seen as an appropriate model for the standard open Weyl chamber of type BCn- The 
cotangent bundle T*c is trivial, whence the phase space of the RSvD system can be identified 
with the product manifold 


= c X = {(A, 0) I A e c, 0 e R”}, 


( 1 . 2 ) 
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that we endow with the symplectic form 


= 2 d6*c A dAc- (1.3) 

C=1 


We mention in passing that the unusual numerical factor in is inserted purely for consistency 
with our earlier works |2Sl ES]- As for the dynamics, it is governed by the Hamiltonian 
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where /x, v and k are arbitrary real parameters satisfying /i < 0 < ix. Also, on these so-called 
coupling constants in this paper we impose the condition uk > 0. As can be seen in [27], this 
additional requirement ensures that the particle system possesses only scattering trajectories. 

Having dehned the models of our interest, now we wish to outline the content of the rest 
of the paper. To keep our present work essentially self-contained, in Section [2] we briefly skim 
through the necessary Lie theoretic machinery and the symplectic reduction background, that 
provide the building blocks of the latter developments. Also, this section allows us to £x the 
notations. Starting with Section [3] we present our new results. Section [3] is the longest and the 
most technical part of our paper, in which we study of the r-matrix structure of the rational 
Cn RSvD model corresponding to the special choice k = 0. Sticking to the Marsden-Weinstein 
reduction approach, in Subsection 13.11 we construct local extensions for the Lax matrix of the 
rational RSvD model. Making use of these local sections, in Subsection 13.21 a series of short 
Propositions and Lemmas allows us to construct a classical r-matrix structure for the C^-type 
model. In this respect our main result is Theorem [101 in which we formulate the r-matrix 
structure in a convenient quadratic form. The resulting quadratic r-matrices turn out to be 
fully dynamical, depending on all variables of the phase space V^. Subsequently, by switching 
to a purely algebraic approach, in Section H] we generalize Theorem [10] to the rational BC^ 
RSvD system with three independent coupling constants. The quadratic r-matrix structure 
of the -BC„-type system is summarized in Theorem [TT] To make this important result more 
transparent, in Theorem [12] we describe the r-matrix structure in a more convenient choice of 
gauge. In this gauge we also provide a Lax representation of the dynamics, as formulated in 
Theorem [T31 Finally, in Section [5] we offer a short discussion on our results and also point out 
some open problems related to the RSvD systems. 


2 Preliminaries 

In this section we overview those Lie theoretic notions and results that underlie the geometric 
construction of the classical r-matrix structure for the rational Cn RSvD system. Our approach 
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is based on the symplectic reduction derivation of the RSvD models, that we also briefly outline. 
In Subsection 12.11 we closely follow the conventions of the standard reference [2S], whereas in 
Subsection l2.2l we employ the notations introduced in our earlier work [23] on the RSvD systems. 


2.1 Lie theoretic background 

Take a positive integer n G N and keep it hxed. Let N = 2n and introduce the sets 


Mn = n} and Nn = {I, ■ ■ ■, N}. 

With the aid of the N x N matrix 

^ _ On In 

^ “ 1 0 

■^n '^n 

we dehne the non-compact real reductive matrix Lie group 

G = t/(n, n) = {ye GL{N, C) | y*Cy = C}, 
that we equip with the Cartan involution 

y^iy-^r. 


Its hxed-point set 

K = {y eG\ 0{y) = y} 

is a maximal compact subgroup of G, having the identihcation K = U{n) x U{n). 
On the Lie algebra 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 


g = u(ti,t!) = {Y 6 sI{N,C)\Y’C+CY = 0} (2.6) 

the corresponding involution 

Y^-Y* (2.7) 

naturally induces the Cartan decomposition 

0 = f©p (2.8) 

with the Lie subalgebra and the complementary subspace 

t = ker('d — Idg) and p = ker('d-|- Idg), (2.9) 

respectively. That is, each element Y E g can be decomposed as 

Y = Y+ + Y_ (2.10) 

with unique components Y^ E t and YL G p. Notice that the Z 2 -gradation fl2.8p of g is actually 
orthogonal with respect to the non-degenerate Ad-invariant symmetric bilinear form 

(,):gxg^M, (W, Rs) ^ MTO). (2.11) 
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To make our presentation simpler, for all /c,/ G Nat we introduce the standard elementary 
matrix G 0l(-^, C) with entries 

{Gk,l)k',l' = {k\l'E'H n)- ( 2 . 12 ) 

Also, with each A = (Ai,..., A„) G we associate the N x N diagonal matrix 

A(A) = diag(Ai,..., A„, — Ai,.. ., — A„) G p. (2.13) 


The set of diagonal matrices 

a = {A(A) I A G 

forms a maximal Abelian subspace in p. Note that in a the family of matrices 


^n+c,n+c) (c ^ f^n) 

forms an orthonormal basis, i.e. {D~, D^) = 6c,d for all c, d G Nn. 

The centralizer of the Lie algebra a inside K is the Abelian Lie group 

M = Zk.(o) = | 6 


(2.14) 


(2,15) 


(2.16) 


with Lie algebra 


m = {diag(ixi,..., iy^, iyi, • • •, iXn) I Xi, • • •, Xn e 
In this Abelian Lie algebra the set of matrices 

(^c,c “ 1 “ 6?i+c,nH-c) (c G 


(2.17) 


(2.18) 


forms a basis obeying the orthogonality relations = —6c,d (c, d G M„). 

Let m-‘“ and a-*" denote the sets of the off-diagonal elements of i and p, respectively. With 
these subspaces can write the refined orthogonal decomposition 


0 =m©m ©a©a. 

In other words, each element Y E g can be uniquely decomposed as 

Y = Y^ + Y^^ + Y, + y;x. 


(2.19) 


( 2 . 20 ) 


where each component belongs to the subspace indicated by the subscript. In order to provide 
convenient bases in the subspaces m.-*- and a"*", for each c G we introduce the linear functional 

Ec- M” ^ M, A = (Ai,..., Aji) I—)■ Ac. (2.21) 

Let us observe that the set of functionals 


TZ+ = {sa ±£fe|l<a< 6 <n}U {2ec | c G 


( 2 . 22 ) 
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can be seen as a realization of a set of positive roots of type Now, associated with the 
positive root 2ec (c G Nn), we dehne the matrices 


^2ec ^^(^c,n+c i ^^n+c,c)- 


(2.23) 


In association with the other positive roots, for all 1 < a < 6 < n we dehne the following 
matrices with pnrely real entries: 


2 T ('b,a i ('n+a,n+b ('n+b,n+a')j 

-±, 1 - 


^£a+£b 2 ^b,n+a i (in+a,b T ^n+b,a) ^ 

together with the following ones with pnrely imaginary entries: 


(2.24) 


2 ^ ^b,a i ('n+a,n+b 4“ ('n+b,n+a} y 

^±,i 


^£a+Sb 2 ^b,n+a i ^n+a,b i ^n+b,a') ■ 


(2.25) 


The point is that the set of vectors forms a basis in the snbspace nx''', whereas the family 

provides a basis in o-*-. Moreover, they obey the orthogonality relations 


(X+^xy) = -<5„,„,<5,,,, and (X’^ X"/) = 


(2.26) 


Note that the family of vectors 

{<’/} = {Df} U {Ay--} (2.27) 

forms a basis in the real Lie algebra u(?7,,n). We mention in passing that it is a basis in the 
complexihcation qI{N,C) = u{n,n)'^, too. 

Next we tnrn to the linear operator 


adA(A): 0 0! ^[A(A), X], (2.28) 

dehned for each A G M”. The real convenience of the basis fl2.27p stems from the commutation 
relations 

adA{A)(-D^) = 0 and adA(A)(A^’') = a(A)Xj’", (2.29) 

where c G N^, a G 71+ and e G {r, i}. Notice that the subspace m-*- © a-*“ is invariant under the 
linear operator adA{A )5 whence the restriction 

adA{A) = adA{A)|m-Lea-L e 0[(m-^ © a-^) (2.30) 

is well-dehned for all A G M"', with spectrum 

Spec(adA(A)) = {±«(A) | a G 71+}. (2-31) 
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The regular part of a is defined by the subset 

Ureg = {A(A) I A e M" and adA(A) is invertible}, (2.32) 

in which the standard Weyl chamber {A(A) | A G c} is an appropriate connected component. 
Note that this Weyl chamber can be naturally identihed with the conhguration space c fll.l|) of 
the rational BCn RSvD system. 

Having set up the algebraic stage, now we turn to some geometric results that are specific to 
the symplectic reduction derivation of the rational RSvD models. First, recall that the regular 
part of p fl2.9l] defined by 


Preg = {kA{X)k~^ I A G c and k G K} (2.33) 

is a dense and open subset of p. It is an important fact that with the smooth free right M-action 

M X {c X K) 3 (m, (A, k)) i-G- (A, km) E c x K (2.34) 


the map 


TT : C X iF -» preg, (A, k) I— )■ kA{X)k ^ 
is a smooth principal M-bundle, providing the identification 

preg = (c X K)/M = C X (K/M). 


(2.35) 

(2.36) 


In the geometric construction of the dynamical r-matrix for the rational Cn RSvD model we 
shall utilize certain local sections of tt with the characteristic properties below. 

Proposition 1. Take an arbitrary point A® G c and let A^°^ = A(A^°^). Then there is a smooth 
local section 

Preg 3Y ^ (e(F), (t{Y)) EcxK (2.37) 

of 71 Id. ,9,5)) . defined on some open subset preg C preg, such that 

A(°) G Preg, (e(A(°)),a(A(°))) = (A^^^, 1), ran(a,^(o)) C m^. (2.38) 

Moreover, under these conditions, at the point A*'°^ the action of the derivatives of e and a on 
the tangent vector SY G p = T^(o)preg takes the form 

e,A(0) • • •, SYn,n) e M" = T^io) c, (2.39) 

e c t ^ TliF. (2.40) 

Proof. Notice that the point (A)°\ 1) G c x iF projects onto A^°^, that is, 

7r(A(°^ 1) = A(°\ (2.41) 

Differentiating vr fl2.35p at (A^°\ 1), let ns observe that for each tangent vector 


(5A © hfc G M” © t = T^(o)C® TiiF = T(;,(o)_i)(c x K) 


(2.42) 
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we can write 


= |Avr(A(o) = A(5A) - [A^^\6k]. (2.43) 

Utilizing the linear operator fl2.3np . it is clear that 

7r*(A(o),i)(^-^ ® = A{6X) - cidj^(o){{Sk)^±), (2.44) 

from where we conclnde that 

ker(7r^(;,(o)^i)) = {0} © m. (2.45) 

Since the snbspace M" ©m-‘“ is a complementary snbspace of ker(7r^(;^(o) i)) in the tangent space 
T(;^(o)_i)(cx A), it is evident that there exists a local section (e, a) 02.371) satisfying the conditions 
imposed in 02.38p . Moreover, by differentiating the equation 

71 o (e, a) = Idp,,g (2.46) 

at the point the relationship 02.441) entails that for all 6Y G p = T^(o)preg we can write 

A(e,^(o) (^U)) - ^^(0) (a,^(o) (<5y)) = (2.47) 

Projecting this equation onto the subspaces a and a-*-, respectively, the formulae for the deriva¬ 
tives displayed in 02.39P and O2.40p follow at once. □ 

To proceed further, we introduce the set of complex column vectors 

= {U G I CU + U = 0 and V*V = A}, (2.48) 

that can be naturally identihed with a sphere of real dimension 2?7, — 1. At each point U G S' 
the tangent space to S can be identihed with the real subspace of the complex column vectors 

TyS = {6V G I C5V + (5U = 0 and {6VyV + V*5V = 0}, (2.49) 

that we endow with the inner product 

(5U, 5v)TyS = Re(((5U)*5n) {SV, 6v G TyS). (2.50) 

Next, we introduce the distinguished column vector E E S with components 

Ea = l and E^+a = -1 (a G M„). (2.51) 

Also, with each vector U G S' we associate the N x N matrix 

e(U) =i/i(UU*-l)+i(/i-z/)CGU (2.52) 

Since the iP-action on S dehned by the smooth map 


K X S 3 {k,V) ^ kV e S 


( 2 . 53 ) 














is transitive, and since k^{V)k ^ = ^{kV) for all k E K and G S', it is clear that the adjoint 
orbit of K passing through the element ^{E) G t has the form 

O = {ki{E)k-^ \keK} = {^(1/) I 1^ G S}. (2.54) 

As is known, the orbit O can be seen as an embedded submanifold of t, and for its tangent 
spaces we have the identihcations 

TpO = {[X,p\\X Ei} {pEO). (2.55) 

In our earlier papers [221 EDI ESI ED] we have seen many times that this non-trivial minimal 
adjoint orbit plays a distinguished role in the symplectic reduction derivation of both the CMS 
and the RSvD systems. In this paper, throughout the construction of a dynamical r-matrix for 
the rational Cn RSvD system, we will also exploit that with the free I7(l)-action 

U{1) X S3 V) ^ E S (2.56) 


the map 


is a smooth principal I7(l)-bundle, providing the identihcation O = S/U{1). 
it is clear that the derivative of ^ takes the form 


(2.57) 
Recalling fl2.52p . 


e*y(5R) = ip{{6V)V* + V{6vy) E T^iyP {V E S, 6V E TyS), (2.58) 


whence it follows that 

ker(e*v) = MiR and (ker(e*y))^ = G TyS \ {6vyV = V*6v}. (2.59) 

Let us also note that for all X G t and R G S' we have XR G TyS and 


e,v(XR) = [X, e(R)] G T^(v)(9. (2.60) 

The last two equations entail that for each 6V E TyS one can hnd a Lie algebra element X G t 
and a real number f G M such that 


6V = XV + tiV. (2.61) 

Having determined the derivative of now we shall work out the derivatives of certain local 
sections, that hnd applications it the latter developments. 

Proposition 2. Let E S be an arbitrary point and define p® = ■^(12^°^) G O. Take a 
smooth local section 

W:d^S, p^W{p) (2.62) 

of ^ (|^.57| ), defined on some open subset O TO, satisfying the conditions 

G a, IR(p(°)) = ran(IR,po)) C (ker(e,v(o)))^. (2.63) 

Then for the derivative of W at the point we have 

IR,po)([X,p(°)]) = XV(°) - ^ (2.64) 
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Proof. It is evident that there is a smooth local section W of the principal f/(l)-bnndle ^ that 
satishes the conditions displayed in fl2.63p . Take an arbitrary tangent vector [X,e Tp(0)O 
generated by some X G and introdnce the shorthand notation 

5hT = lV,,(o)([X,p(°)]) e (ker(e*v(o)))^. (2.65) 

By taking the derivative of the relationship ^ oW = Id^ at the point we find that 

e*v(0)(5W^) = e*v(0) o lX,po)([X,pW]) = [X,p(0)] = e*v(0)(^VW), (2.66) 

therefore 5W — XV^°^ G ker(,^^y(o)). However, due to fl2.59p we can write that 

5W = XV(°) + (2.67) 

with a unique real number x. Its value can determined by the fact that the tangent vector 5W 
belongs to subspace (ker(,^^y(o)))■*■, leading to the formula 02.641) . □ 

2.2 The rational RSvD model from symplectic reduction 

Based on our earlier results, in this subsection we review the symplectic reduction derivation of 
the rational RSvD system. The surrounding ideas and the proofs can be found in |2S]. An 
important ingredient of the symplectic reduction derivation of the RSvD system of our interest 
is the cotangent bundle T*G of the Lie group G 02.3p . For convenience, we trivialize T*G by the 
left translations. Moreover, by identifying the dual space g* with the Lie algebra g fl2.6p via the 
bilinear form fl2.1ip . it is clear that the product manifold V = G x g provides an appropriate 
model for T*G. For the tangent spaces of the manifold V we have the natural identifications 

~ ® ^yQ — TyG © g ((p, Y) G P), (2.68) 

and for the canonical symplectic form u G PL^{V) we can write 

© AX, 5y © 5Y) = (p“^Ap, 5Y) - {y~^5y, AY) + {[y~^Ay, y~^6y], Y), (2.69) 

where (p, X) G P is an arbitrary point and Ay © AX, 6y © 6Y G TyG © g are arbitrary tangent 
vectors. An equally important building block in the geometric picture underlying reduction 
derivation of the RSvD model is the adjoint orbit O fl2.54p . Of course, it carries the Kirillov- 
Kostant-Souriau symplectic form G that can be written as 

©(|A', p\, [Z, p\) = (p, [A, Z]) (p 6 O, A, Z 6 t). (2.70) 

Making use of the bundle fl2.57p and the equations fl2.60p and fl2.6ip . one can easily see that 

= 2p\m{{8Vy8v) (X G 8V, 8v G TyS). (2.71) 

Now, by taking the symplectic product of the symplectic manifolds {V,uj) and (C>,a;®), we 
introduce the extended phase space 

(pext,^ext) _ (P X C>,a; + a;^). (2.72) 
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To describe the Poisson bracket on this space, for each smooth function F G at each 

point u = (y, F, p) G we dehne the gradients 

V^F(m) g g, V®F(m) g g, V^F(m) g TpO (2.73) 

by the natural requirement 

FU5y®5Y®[X,p]) = (y^F{u),y-^5y) + (y^F{u),5Y) + (V^F(n),X), (2.74) 

where 6y G T^G, 6Y G g and X G 6 are arbitrary elements. Now, one can easily verify that the 
Poisson bracket on induced by the symplectic form can be cast into the form 

{F, HY-\u) = (V^F(«), Vm{u)) - (V3F(n), V^H{u)) 

-{[V^F{u),V^H{u)lY)+u^{y^F{u),V^H{u)), 

for all F,He G°°(P®^*). To proceed further, let us note that the smooth map 

: (X X X) X ^ {{kL, kR), (y, Y, p)) ^ {kLyk]^, knYk^, kLpkY) (2.76) 

is a symplectic left action of the product Lie group K x K on the extended phase space 
and it admits a K x iP-equivariant momentum map 

jext. pext ^ ^ Y, p) ^ YyYy-^)+ + p) © (-F+). (2.77) 


As we proved in [25], the rational G„ RSvD model can be derived by reducing the symplectic 
manifold P®’^* at the zero value of the momentum map 

Let us recall that the standard Marsden-Weinstein reduction consists of two major steps. 
At the outset, we need control over the level set 

£0 = (P"*)-'({0}) = {m G P""* I r^\u) = 0}, (2.78) 


that turns out to be an embedded submanifold of P®^* fl2.72p . However, to get a hner picture, 
we still need some more background material. First, for each a G we dehne the rational 
function 


c 3 A I —y Xq(A) 



n 

d=l 

{d^a) 




2i/i A 
Aa + Ad / 


G C. 


(2.79) 


Also, we need the vector-valued function F : P^ —)■ with components 

Fa = (^^"\ZaY and Fn+a = (Y^Fa\Za\~^ (oGN^), (2.80) 

that allows us to introduce the function A: P^ —)■ exp(p) with the matrix entries 


2ipFaFb 


Aa^b 
Aa,n+b An+b,a 


A: 


2i/i + Aq — Af, 

2ipFr^Fn+b 


2ipFn+aF n+b 
n+a,n+b — W: 7 ; ^) 

2l/i - Aa + Afe 
i(/i - z/) 


+ 


2i/i + Afl + Afo i/x + A 




(2.81) 
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where a, 6 G As we have seen in [25], function A provides a Lax matrix for the rational Cn 
RSvD model with the two independent parameters fi and u. Next, let us consider the smooth 
function V; —)■ S defined by the equation 

V = (2.82) 

and also introduce the product manifold 

= V^ X {K X K)/U{1),, (2.83) 

where 17(1)* stands for the diagonal embedding of 17(1) in the product group K x K. Having 

equipped with the above objects, now we are in a position to provide a convenient parametriza- 
tion of the level set £0 fl2.78p . Indeed, in [25] we proved that the map 

(2.84) 

defined by the assignment 

{X,9,{vl,Vr)U{1)*) ^ {riLA{X,9)^ri]^^,riRA{\)ri^\r]L^{V{X,9))ril^) (2.85) 

is a smooth injective immersion with image = Hq- Moreover, in [25] we also proved 

that gives rise to a diffeomorphism from onto the embedded submanifold ilo- In other 
words, the pair T^) provides a model for the level set £0 112.78p . 

To complete the Marsden-Weinstein reduction, notice that the (residual) K x 77-action on 
the model space A4^ 112. 83 p takes the form 

(kL, kn) ■ (A, 9, {r]L, t]r)U{1)^) = (A, 9, {kiVL, kRT]R)U{l)^), (2.86) 

thus the orbit space Ai^/{K x K) can be naturally identified with the base manifold of the 
trivial principal (77 x 77)/17(l)*-bundle 

^R, {XAAvl,Vr)U{1).) ^ (A,0). (2.87) 

Now, the crux of the matter is the relationship 

= (T^)*a;""\ (2.88) 

that we proved in |2S] by applying a chain of delicate arguments. Therefore, for the symplectic 
quotient in question we obtain the identification 

(P""V/o(i7 X (2.89) 

Finally, note that the 77 x 77-invariant function 

Y ; ^ M, (y, Y, p) ^ Ui{yy*) (2.90) 

survives the reduction, and by applying straightforward algebraic manipulations one can verify 
that the corresponding reduced function coincides with the Hamiltonian of the rational 
RSvD system with two independent coupling parameters p and z/, that can be obtained from 
the RC^-type Hamiltonian lll.4p by setting k = 0. 
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3 Dynamical r-matrix for the (7r^-type model 

Building on the symplectic reduction picture outlined in the previous subsection, our goal is to 
construct a classical r-matrix structure for the C,i-type rational RSvD system with two inde¬ 
pendent coupling parameters. In the context of the CMS models, this geometric approach goes 
back to the work of Avan, Babelon, and Talon |2I] . Eventually, in our paper [32] , we succeeded 
to construct a dynamical r-matrix for the most general hyperbolic BCn Sutherland model with 
three independent coupling constants, too. It is worth mentioning that the surrounding ideas 
proves to be fruitful in the broader context of integrable field theories as well. For a systematic 
review see e.g. [33] . 

As we have seen in [23], the eigenvalues of the Lax matrix A 02.811) do commute, whence 
it follows from general principles that A obeys an r-matrix Poisson bracket (for proof, see e.g. 
[saias]). However, we wish to make this r-matrix structure as explicit as possible. For this 
reason. Subsection [Q is devoted to the study of certain local extensions for the Lax matrix 
of the rational Cn RSvD model. As it turns out, these local extensions are at the heart of the 
construction of the dynamical r-matrix structure for the RSvD system, that we elaborate in 
Subsection 13.21 

3.1 Local extensions of the Lax matrix A 

The backbone of our reduction approach is the construction of the so-called local extensions of 
the Lax operator A fl2.8ip . that we wish to describe below. For this reason, take an arbitrary 
point 

e (3.1) 

and keep it fixed. Clearly the point 

g(0) _ ^^(0)^ ^(0)^ l)f/(l)*) e (3.2) 

is one of the representatives of (A^°\ 6*^°^) in fl2.83p . that is, 7r'^(s®) = (A^^^ 6*®). Moreover, 
let us introduce the shorthand notations 

= M(A(°), = J^(A(°), = V(A(°), 0(°)), (3.3) 

together with 

2/(o) = (MW)T rW = = A(A®), pW=e(V(°)). (3.4) 

Corresponding to 5*-°^ fl3.2p . in the extended phase space we also introduce the reference point 

^(0) _ xK(s(o)) _ (^(0)^ r(o), pW) G (3.5) 

Now, associated with the elements given in fl3.31l3.4p . let us choose a local section (e, a) of vr 
02.351) . and also a local section IF of ^ 02.571) . as described in Propositions [Hand [21 respectively. 
Upon defining the open subset 


fl = {R G g I R_ e preg} R g, 


(3.6) 
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it is clear that 


^: G X 0 X (!) ^ C^, (y, F, p) ^ a{Y_)-^y*W{p) (3.7) 

is a well-defined smooth function. Due to the conditions imposed in the equations 02.381) and 
02.631) . at the point 03. 5 p for the first n components of d' we have 

{aeNn). (3.8) 

Since these components are strictly positive, there is an open subset C GxgxO containing 
the distinguished point such that for all a G the map 

^ (u) 

m,:P®^‘^f/(l), (3-9) 

is well-defined and smooth. Let us keep in mind that by construction ma{u^^^) = 1. 

Now we are in a position to define those group-valued functions that play the most important 
role in the construction of a dynamical r-matrix for the rational G„ RSvD system. First, making 


use of the functions rria fl3.9p . we build up the M-valued function 

M I—)■ diag(mi(M),...,m„(n),mi(M),...,mn(M)), (3.10) 

which satisfies = 1. Next, we introduce the iF-valued functions 

(p,F,p)^a(F_), (3.11) 

(p: P®"* k{u)m{u). (3.12) 

Clearly we have = 1. Finally, we define the functions 

A(o):pext^G, iy,Y,p)^y*y, (3.13) 

A:P®"*^G, u^p{u)-^A^°\u)p{u). (3.14) 


Notice that at the point fl3.5p we have A^°^(m^°^) = A(n^°^) = Having equipped with 
the above objects, now we can formulate the central result of this subsection. 

Lemma 3. The G-valued smooth function A l\3.14\ l ® local extension of the Lax matrix A 

H2.81\) around the point in the sense that AfvJ^) = A(A®,6*(°)) and 

(3.15) 

Proof. It is enough to verify fl3.15p . For, take an arbitrary point 

(A,9,(i7i,i7B)C/(l),) e (T«)-‘CP“) (3.16) 

with some A G c, 0 G M" and pl, pr G K. Also, for brevity we define 

t< = (9,F,p) = T«(s)eP“‘. (3.17) 
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Recalling fl2.35l) and (12.851) . it is clear that 

V = rjRA(\)T/^‘ = ir(A, ijr) £ p„e- 


(3.18) 


On the other hand, since V G g, we have YL G preg- Thus, making use of the local section (e, a) 
introduced in fl2.37p . we see that 


Y = Y_ = 7r(e(y_),cr(yi)) = 7i{e{Y_), k{u)) 


(3.19) 


also holds. Recalling fl2.34p . the comparison of fl3.18|) and fl3.19p yields that there is a unique 
element 

fh = diag(mi,..., m„, fhi ,..., m„) G M (3.20) 

such that 

(A,r?K) = (e(Y'_),fc(M)m). (3.21) 

Next, remembering the parametrization fl2.85p . we can write 

aViViX, 9)) = 0))hZ' = P e a. (3.22) 

However, utilizing the local section W introduced in fl2.62p . we also have i{W{p)) = p, whence 
by (12.561) we can write that 

PlV(A, P) = (p) (3.23) 

with some constant G M. From the above observations it readily follows that 

e'’/’VF(p) = plA{X, O) = {vT^VrHK O) = {y*)-"k{u)mF{X, 9), (3.24) 

from where we get = rhJ^{X, 9). Componentwise, for each a G we can write 

e^Ha{u)=maJ^a{X,9), (3.25) 

thus the relationship |\['a(M)| = J^a{X,9) and also 




MX,9) |T„(m)| ^ 


(3.26) 


are evident. In other words, fh = e^'^m{u), whence from (13.211) and (I3.12p we conclude that 

rjji = k{u)m = k{u)m{u) = e^'^(p{u). (3.27) 

Now, turning to the functions (I3.13|) and (I3.14|) . notice that 

(m) = y*y = r]RA{X, 9)M, (3.28) 

so from (13.271) we infer that 

A o T^(s) = A{u) = ip{u)~^ A^^\u)(p{y) = Al(A, 9) = Ao 7r'^(s). (3.29) 

Since s (I3.16P is an arbitrary element of (T^)“^('P®’^*), the Lemma follows. □ 
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3.2 Computing the r-matrix 

The natural idea impregnated by Lemma [3] is that the Poisson brackets of the components of 
the Lax matrix A fl2.81jl can be computed by inspecting the Poisson brackets of the components 
of the locally defined function A fl3.14jl . Indeed, since we reduce the symplectic manifold 
fl2.72l) at the zero value of the K x iP-equivariant momentum map fl2.77p . and since the local 
extension A is (locally) K x iP-invariant on the level set £o ll2.78p . using the St. Petersburg 
tensorial notation we can simply write 

{A ® = {A ® (3.30) 

However, for the function fl3.13p we clearly have 

{^(0) ® ^(0)}ext _ 0, (3.31) 

that is, obeys an r-matrix bracket with the trivial zero r-matrix. Therefore, due to the 
relationship A = if fl3.14p . it is clear that A also obeys a linear r-matrix bracket 

{A f A}®"* = [7i2, a ® 1] - [rai, 1 <8) A] (3.32) 

with the transformed r-matrix 

7 i 2 = ¥^rV2 ^ A®}®^vr^ + 7 ^ 2 }“VrV^\ Af (^1(^2. (3.33) 

Now, recalling that = 1, from the relationships fl3.3np and fl3.33p we infer that for the 

Lax matrix A we can write 

{Ai,A 2}''(AW,0(°)) = [ri2(A(°),0(°)),Ai(A(°),0(°))] - [r2i(A(°),0(°)),A2(A(°),0W)] 
with the r-matrix 

ri2(A<"',»<“') = »<“>)]■ 

However, since ip = km fl3.12p . Leibniz rule yields 

together with 

{ipup>2r\u^^^) ={kuk2r\u^°'>) + {k,,m2r\u^°'>) 

Thus, in order to provide an explicit formula for the above r-matrix 03.351) . we still have to work 
out the Poisson brackets appearing in 03.361) and 03.371) . However, recalling 02.75p . it essentially 
boils down to the computation of the gradients 02.731) of the components of the functions A^°\ 
k and m. We accomplish these tasks in the following series of short Propositions. 


(3.34) 

(3.35) 

(3.36) 

(3.37) 
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(3.38) 


Proposition 4. At the reference point \3. for all matrix v G Qi{N, C) we have 


(V^Re(tr(vA(°))))(«(°)) = 
(V^Im(tr(uA(°))))(M®) = 


2 

2 i 


whereas the remaining gradients o/Re(tr(i;A®)) and Im(tr(i;A^°^)) are trivial, i.e. 

(V0Re(tr(nA(°))))(n(°)) = (V0Im(tr(nA(°))))(n(°)) = 0, 

(V^Re(tr(i;A*'°^)))(n‘'°^) = (V^Im(tr(i;A*'°^)))(n‘'°^) = 0. 

Proof. Take an arbitrary tangent vector 

^u = 8y®6Y® [X, G T,(o)P“* 


(3.39) 

(3.40) 

(3.41) 

(3.42) 


with some Lie algebra element X G t. By neglecting the second and the higher order terms in 
the small real parameter t, one can easily find that 


A(°)(-u(o) + tAu + • ■ ■) 

= .41”' + t (.4<“>((!/'“>)-'i9) + + ... , 

from where we infer that 

Re(tr(nA^°^))^„(o)(An) = |^Re(tr(nA^°^(n^°^ + tAn + ■ ■ ■))) 


(3.43) 


t=o 


= tr 


(n + ^;*)A(°) + (A(°))-H^^ + ^^*) 


(d®) ^^ 2 /) • 


(3.44) 


Similarly, one obtains immediately that 

Im(tr(„A'»))),„<.,(Au) = tr |' (f " - v>) ^ (3 45 ) 

Since the N x N matrices appearing on the right hand side of both 03.381) and 03.391) do belong 
to the Lie algebra 0 02.61) . by the definition of the gradients 02.74p the Proposition follows. □ 

Proposition 5. At the distinguished point / Id. ,51) . for all Lie algebra element v E Q we have 
(V®Re(tr(nfc)))(M^°^) =-^(ad^co))”^ ((n — n*)^x) and (V®Im(tr(n/c)))(n*^°^) = 0, (3.46) 

while the remaining gradients of R.e{tr{vk)) and Im(tr(nfc)) are trivial, i.e. 

(V‘^Re(tr(nfc)))(M‘'°^) = (V'^Im(tr(nfc)))(M‘'°^) = 0, 

(V‘^Re(tr(nfc)))(M*^°^) = (V^Im(tr(n/c)))(M*^°^) = 0. 
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(3.47) 

(3.48) 















Proof. As in the proof of Proposition 01 take an arbitrary tangent vector Am as given in fl3.42p . 
By applying a first order expansion on k fid.lip in the small real parameter t, Proposition [T] 
tells us that 


+ tAu + •■■) = cr((A*-°^ + t6Y)_) = (j(A^°^ + t{6Y)^) = 1 + t6a + • • • , 

where ^ 

Sa = -(ad^(o))“^((5F)„±) G m-^. 

Since for all m G 0 we have 

tr^vk) + tY{v*k~^ 


Re{tT{vk)) = 
it is now evident that 


and Im(tr(Mfc)) = 


tT{vk) — tT{v*k 




2 i 


+ iA« + ...))) = + itr 


V — V 


5a + ■ ■■ . 


(3.49) 

(3.50) 

(3.51) 

(3.52) 


Recalling the Cartan involution fl2.7p we can write v — v* = v + i9(v) G therefore 


(Re(tr(MA;))),„(o)(AM) = (n +'^(m), (ad^m) ^((5R)o 


(3.53) 


= -(^(ad^(o)) 

In complete analogy with fl3.52p . for the function Im(tr(M/c)) we can write the expansion 


Im(tr(MA:(M^°^ + tAu + ■••))) = 


tr(M) — tr(M*) 
2 i 


+ ftr 


v + v* 

2 i 


5a + ■ ■■. 


(3.54) 


However, since v + v* = v — i9(v) G p, and since the subspaces t and p 02.911 are orthogonal with 
respect to the bilinear form 02 . 1111 . we conclude that 


(Im(tr(Mfc))),„(o) (Am) = — (m - ■^(m), 5a) = 0, 


thus by recalling 02.7411 the Proposition follows. 


(3.55) 

□ 


To hnd the gradients of the components of m O3.10p we find it convenient to introduce the 
auxiliary function 


r: C\{0}^17(1), 

It is clearly smooth, and at each point x > 0 for its derivative we have 


= 1 - 


,Im(t(;) 


X 


{weC^T,{C\{0})). 


Note that with the aid of r the function rria 03.911 can be simply written as 

ma = To Ta. 


(3.56) 


(3.57) 


(3.58) 
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Also, for each a G we introduce the N x N matrix 

^ _.e,(^(o))* + C(^W)e:C 

Ca 1 2 

where G C'^ denotes the column vector with components 

(cfc)/ = 5k,I {k, I G Nat). 


(3.59) 


(3.60) 


As a matter of fact, the above matrix (a belongs to the Lie algebra g. Moreover, utilizing the 
basis fl2.27p . we can write 

2C. =V2Ti<‘>D: - V2Re(^nl)(Xii + 

a —1 

+ E - xgij - Re(4”+c)(v;p. + xziJ 

C=1 

(3,61) 

n 

+ E (^i”’(M'5. + A'EeJ - R<!(4»i)(AE.. + x;aj 

c=a+l 

-Im(JF»i)(E’k + .V:P.)). 

Proposition 6. Take an arbitrary a G N^, then at the reference point Lt. ,51) the gradients 
1^2. 731 ) of the function Re(ma) are all zeros. However, for the imaginary part of the function 
li3.9\) we have the non-trivial formulae 


(V^(Im(mJ))(M(°)) = ^Ca e g. 


ra 

(V®(Im(mJ))(M(°)) = ^(ad^(o))“^((Ca)mA) e a^. 


T 

a 


(V^(Im(m,)))(M^^0 = e.v(o)(^K) e Txo)0, 


,(o)^ = 


where 


,,, C(MW)ie,-(M(°))^e„ , V® ^ 

5 Va = - m -^ ^ 7 v( 0 )h. 




(0) 


2piN 


(3.62) 

(3.63) 

(3.64) 

(3.65) 


Proof. First, take an arbitrary tangent vector Sy G Ty{o)G. Recalling fl3.7p . and the conditions 
fl2.38p . fl2.63p . for small values of the real parameter t we can write the hrst order expansion 


(3.66) 


= a((E(°))_)-^(2/(°) + % + ••• )*hF(p(°)) = - tC{y^^^)-\5y)CH^^^ + • • • , 

thus for each a G N„ we have 

+% + ••■, p(°)) = +twa + --- (3.67) 
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with Wa = —tr(C(|/*^°)) ^{5y)CF^'^'>e*^) G C. Recalling fl3.57p . we can write 
+ %H- +tWa^ -) = 1 + 


+ • • • , 


(3.68) 


thus clearly (V'^(Re(ma)))(-u‘^°^) = 0. Moreover, by inspecting Wa and fl3.59p . we obtain 

Ini(M;a) tr(Ca(2/^°^)“^(52/) 


(Ini(ma))^„(o) ((5|/ © 0 © 0) = 


T' 

-J C 


( 0 ) 


-pv 

-j a 


(0) 


(3.69) 


so (13.6211 also follows immediately. 

Second, take an arbitrary tangent vector SY G g = Ty(o)g. According to Proposition [H for 
small values of f G M we have the hrst order expansion 

a((F(°) + t6Y)_) = a(A(°) + t{6Y)_) = l + t6a + -- - , (3.70) 

with the Lie algebra element 6a G m-*- displayed in (13.50^ . Therefore, 

^(^( 0 ), y( 0 ) ^ p(°)) = cr((y(°) + t6Y)_)-\y^^^yW{p^^^) = - t{6a)T^^'> + • • • , (3.71) 

and so for each a G we can write 

T,(|/W, + t5Y, p(o)) = +tw'^ + ... (3.72) 

with = — ti{{6a) el) G C. Utilizing r fl3.56p and its derivative fl3.57p . we obtain 

ma{y^^\ + t6Y, = r(J'® +twl + ---) = 1 + + • • • , (3.73) 

J~a 

therefore (V®(Re(ma)))(-u(°^) = 0 is immediate. Remembering fl3.59p it is also clear that 

= (Ca,M = - ^(Ca)mi, (adA(0))"^((<^^)a^)) = ( (ad^(o) )"^ ((Ca)mi ), • (3-74) 

Thus, by combining fl3.73p and (13.7411 . we end up with the formula 

Im«) 


(lm(ma))*„(o)(0 © © 0) = 


F 

■J c 


( 0 ) 


Tc 


( 0 ) 


(ad^(o))-^((Ca),n^),5R), (3.75) 


that readily implies (I3.63p . 

Third, take an arbitrary X G L Remembering (I2.38p . (I2.63p . and (13.7p . notice that 

+ f[X,p(°)] + ■■■) 

whence for all a G we can write 

+ f[x, pW] +...) = + • • • 


(3.76) 


(3.77) 
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with the complex number 


m;" = 




( 0 ) 




-tr(XV(°)(V(°))*). 


It readily follows that 


m, 


+ t[X, p(°)] + •••) = + to" + ...) = 1 + ti 


, Im(ta 


to 

a 


( 0 ) 


+ 


from where we get at once that (V‘^(Re(ma)))(M^°^) = 0 and 


(Im(ma))*„(o) (0 © 0 © [X, p(°)]) = 


Im«) 


TTl' 
•J a 


(0) 


(3.78) 


(3.79) 


(3.80) 


At this point notice that iV*'°^(V*'‘^^)* G t. Therefore, recalling fl3.59p . we can write 


77 ( 0 ) 

Im«) = -tr((^(°))Ka(^®)“ 2 X) + ^tr(iV(°nV^°^)*X) 


— 77 ( 0 ) 

^ n 


T' 

■J c 


( 0 ) 


/ i iV 


(3.81) 


Now, one can verify that the column vector 514 displayed in fl3.65|) does belong to the tangent 
space Ty{o)S 02.491) . Furthermore, recalling 02.58P we hnd that 

1 / / . inlN 1 ^ . mix 1 \ i ln\, ^ 

(3.82) 


x; 


( 0 ) 


(.4™)iCa(.4™)-4+ ^V<"'(V‘”>)- = 


□ 


thus the relationship 03.641) also follows. 

Having the necessary gradients at our disposal, now we are ready to work out the tensorial 
Poisson brackets appearing in 03.36P and 03.37p . 

Lemma 7. At the point \3.^) we can write 


{k ® («(“)) = aP [ 


X+’^ © X; 


a 


(A(«)) 


E 

a,e 

n 


X+’^ © X; 


a 


(A(o)) 


•^2 1 


{m ® = -4”> ( 5^ 4 ® S<“> ) - ( E "t ® S't" 14°’, 

where for each a ENn we have 


(3.83) 

(3.84) 


va=l 


^a=l 


= 


x/2Xi' 


( 0 ) 


V2Re(toh)W2 

2A<"' 

a—l / -^{0) -y—,i 
^ / ‘X c 

C=1 

n 

^E 

c=a+l 


Re(41).WE , Im(to+c).Wt.. 


© 0 ) _ © 0 ) 
'Ac -Aa 


A® + A® 


+ 


a(°) + Ai°) 


(3.85) 




Im(41)W4A 1 


Ai“' + A<“' 

A® ©A® )l 
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Proof. To prove f|3.83p . we utilize the family of iV x matrices {u/} defined in fl2.27l) . that 
forms a basis in the complex linear space gl(iV, C). Recalling Proposition 13 we see that among 
the members of the basis {vj} only the vectors (a E 7Z+, e E {±}) generate non-trivial 
gradients of the form 

V(Re(tr(A'+->A^)))(,.'»)) = (3.86) 

Let C g[(iV, C) denote the dual basis of {u/} provided by the trace-pairing of gl(iV, C). 
Recalling Proposition H] and the explicit expression of the gradients fl3.86p . the Poisson bracket 
formula 02.751) allows us to write 




i,j 




a,e J 


E 


Q:,e 


a(A(o)) 


a(A(o)) 


(3.87) 


so 03.831) follows immediately. 

Making use of Propositions S] and El let us notice that the Poisson bracket formula 02.75P 
yields 


n N 

[m f = - ^ + e„+a,n+a) ® ek,l 

a=l k,l=l 
n 

= - ^ T>+ 0 \/2 (M® (V9(Im(m,)))(n(°)) + (V®(Im(m,)))(n(°V^°^) • 

a=l 


(3.88) 


Therefore, by projecting the Lie algebra element C,a 03.6 ip onto the subspace rm'-, the application 
of 03.63P immediately leads to 03.84p . □ 

Lemma 8 . At / (,?.,5]) we have the trivial Poisson bracket {k f /c}®^*(-u^°^) = 0, whereas 




a=l 


where for each a ENn we have 


2 ^( 0 ) _ 


1 


/ y2Re(ji"y)A'i; 




2 a1 “’ 


E 

C=1 

n 

E 

c=a+l 




x(0) 00) 

Ac Aa 




+,i 

Sc-\-£a 




-+,r 

^€c-\-Sa 


A‘”> + a1”> 

,( 0 ) 


Af + a1“' 


Re(X."y).Y.y.. ^ Im(Afi).Y;+.. 


+,r 


00 ) 00 ) 
Xa Xc 


Ai°) + Ai°) 


Ai°) + a(°) 


(3.89) 


(3.90) 
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Proof. Working with the basis {n/} fl2.27p of gl(iV, C), from Proposition [5] we see that 

V®(Re(tr(n//c)))(M®) G a-*- C p. (3.91) 

Keeping in mind the orthogonal Z 2 -gradation fl2.8p . the Poisson bracket formula fl2.75p gives 
rise to the relationship 


{tr (u/Zc), tr (uj/c) ) 

= ([V0(Re(tr(n/fc)))(M(°)),V0(Re(tr(njfc)))(M(°))],A(°A = 0, 


(3.92) 


thus the equation {k f = 0 follows immediately. 

As we have already observed at the beginning of the proof of Lemma [71 among the members 
of the basis {u/} only the vectors v = generate non-trivial gradients for the component 
functions tr(nZc) at the point Utilizing these gradients fl3.86p . Propositions 0 and E] allow 
us to write 

n 

{m ® Zc}®^*(m®) = '^'^{ma,tT{X;f’^k)y\u^^'^){ea,a + en+a,n+a) ® (-A+’") 

a=l Q;,e 


EE (V^(Im(m,))(n(°)), V®(Re(tr(X+’^A;)))(n(°))) V2Df O 

a=l a,e 

/n ~ 

E ® E<C“' 

a=l ^^ 


(3.93) 


+ ,£ 


a,t 


P} ~ 


a=l 


TTV 
^ a 


Remembering the explicit formula of (a 113.611) and the commutation relations fl2.29p . the Poisson 
bracket 03.891) also follows. 

Lemma 9. At the distinguished point n® 


we have the Poisson bracket 
(0) / n+ n+ 


□ 


(3.94) 


with the coefficients 


= 


l<a<b<n 


+ 


\( 0 ) x( 0 ) 

Aa \ 


Proof. Using the antisymmetry of the Poisson bracket, we hnd that 


(3.95) 


{m f m}“*(M(°)) = ^ {rua, mhY^\u^^'^){ea,a + e^+^.n+a) ® {eb,b + en+b,n+b) 

(3.96) 

= -2 {m„m,r*(n(°))(Zl+®Zl+-Zl+®Zl+). 

l<a<6<n 
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To proceed further, let us choose arbitrary a, 6 G satisfying a < b. Notice that the Poisson 
bracket formula fl2.75l) naturally leads to the expression 


+ ((V0(lm(m,)))(u(°)), (V''(Im(mb)))(u(°))) (3.97) 


However, by utilizing Proposition [HI each term on the right hand side of the above equation 
can be cast into a fairly explicit form. Starting with the first term, the application of fl3.6ip 
gives rise to the relationship 


1 [ , Re(4°h)fe(d“7) 

+ Ai°) + A® y 

Keeping in mind that a < b, a similar argument provides 
((V0(Im(m,)))(n(°)),(V''(Im(mb)))(n(°))) 

^ 1 Re(J-lt)Re(J-l%) _ Im(J-(t)Im(J-(;j \ (3.99) 

V Ai°^ - A[°^ Ai°^ + Af Ai°^ + A® J ’ 


Now, let us turn to the third appearing in 03.971) . Utilizing the concise formula 02.7ip for the 
symplectic form 02.70p . the application of the equations 03.641) . 03.65P and 02.8ip yields that 


(Im(ma))) (m(°) ), ((Im(mfe))) (m®)) 


= 2fj.lm{{5VaySVb) = 2plm 




A 


(0) 

a,b 




— '^|(V(0))('C*V(0) (<^^( 1 ), (^^y{0) (514)) 

1 Ai°^ - A^ 

2 (A®-A ®)2 +V' 


(3.100) 


Now, by simply putting together the above equations, the Lemma follows at once. □ 

At this point we are in a position to provide an explicit formula for the r-matrix 03.35p . 
Remembering 03.36p . let us notice that Lemma [7| itself implies that r is in fact linear in A, 
having the form 

ri 2 (A<»',«'“>) = + 4°’(pr2)<“’ (3.101) 

with the g (8) g-valued functions 


(pfj)'"’ 


^ X*' » A--.- 
^ a(AI»)) 


n 

+^D„+®sf ±-{p®pr‘(«<“>). 

a=l 


(3.102) 


Recalling 03.37p . the above expressions can be further expanded. Indeed, by simply plugging 
the formulae displayed in Lemmas [8] and [9] into O3.102p . we may obtain explicit expressions 
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for both r. However, since r is linear in A as dictated by fl3.10ip . the linear r-matrix 

Poisson bracket fl3.34p can be cast into a quadratic form. Also, since the point fl3.ip 

we hxed at the beginning of Subsection 13. II was an arbitrary element of the zero superscripts 
become superfluous and can be safely omitted. With the usual conventions for the symmetric 
and the antisymmetric tensor products, 

XyY = X®Y + Y®X and X ^Y = X®Y -Y ®X, (3.103) 

we end up with the following result. 

Theorem 10. The Lax matrix A ^2.81\] of the rational Cn RSvD model with two independent 


coupling parameters obeys the quadratic Poisson bracket 

{A f = cii2AiA2 + Aibi2A2 — A 2 C 12 A 1 — AiA2di2 (3.104) 

with the g <8) Q-valued dynamical structure coefficients 

Q12 = -^ ^-Da A (Aa + T(j) + ^ TafiDfADf, (3.105) 

a,e ^ '' a=l l<a<b<n 

+ Y, (3.106) 

^ a=l l<a<b<n 

ci2 ^ ^ ’Iff*' - ^ -S'- - A r„) + A r>+, (3.107) 

di2 = ^ + x; A (g. - r.) - 5^ A r>+, (3.108) 

o,e ' a=l l<a<b<n 

where the constituent objects are defined in Lemmas 00 and 0 

Proof. Due to 03.1011) . the Poisson bracket 03.34p takes the quadratic form O3.104p with 

«12 = Pu - pti^ bi2 = -pt2 - PAi Ci2 = -pti - Pu^ di2 = Pu “ Pfv (3.109) 

Remembering 03.1021) . 03.37p . and the explicit formulae displayed in Lemmas [ 8 ] and [9l the 
Theorem follows. □ 


We conclude this section with an important remark. Since the quadratic structure matri¬ 
ces 03.105113.108p are derived from an r-matrix linear in A as described in O3.10ip . from the 
relationships O3.109p it follows immediately that they satisfy the consistency conditions 

0-21 = —0,12, d2l = —di2, ^21 = Ci 2 , 0.12 + ^12 = C 12 -|- <il 2 - (3.110) 

The above observation can be paraphrased as follows. If a Lax matrix A obeys a tensorial 
Poisson bracket 03.341) . and if the governing r-matrix is itself linear in A as in O3.10ip . then 
the tensorial Poisson bracket can be rewritten as a quadratic bracket O3.104p with quadratic 
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structure matrices obeying the consistency conditions fl3.110|) automatically. It is a nice, but 
essentially trivial algebraic fact that the converse of this statement is also true. Indeed, suppose 
that a Lax matrix A obeys a quadratic Poisson bracket fl3.1U4|) with coefficients satisfying 
fl3.1injl . Under these assumptions the quadratic bracket can be cast into a linear form fl3.34jl . 
More precisely, the governing r-matrix can be written in the form of fl3.101jl with 


+ 0,12 + 012 1 _ di2 — bi2 — Ci2 — Ui2 

P12 = - 2 - ^‘2 =- 2 - 


( 3 . 111 ) 


where Uu is an arbitrary g V g-valued function on the phase space, i.e. it obeys the symmetry 
condition U 21 = U 12 . This observation plays a crucial role in the developments of the next 
section. 


4 Classical r-matrix structure of the BCn-type model 


Utilizing a symplectic reduction framework, so far we have studied the classical r-matrix struc¬ 
ture for the rational Cn RSvD model with two independent coupling parameters and u. 
However, to handle the BCn-tjpe model as well, in this section we slightly change our point 
of view. Switching to a purely algebraic approach, we shall generalize Theorem [TU] to cover 
the most general rational BCn RSvD model with three independent coupling constants. As an 
added bonus, at the end of this section we will provide a Lax representation of the dynamics, 
too. 

To describe the Lax matrix of the rational BCn RSvD system with the additional third real 
parameter k, we need the functions 


\/X + \/x2 ^2 1 N . 1 1 

OL[x] = -;=- and p{x) = ik- 


\/X + yfx^ -|- 

where x G (0, cx)). Also, with each A = (Ai,..., A„) G c we associate the group element 


(4.1) 


h{\) 


diag(Q:(Ai),...,Q;(A„)) 

-diag(/3(Ai),...,/3(A„)) 


diag(/3(Ai),... ,/3(A„)) 
diag(Q:(Ai),...,Q:(A„)) 


G G. 


(4.2) 


In 


we proved that the smooth function A: —)■ G defined by the formula 


A(X,e) = h(X)-'^A(X,e)h(X)-' ((A, 9) e P®) 


( 4 . 3 ) 


provides a Lax matrix for the rational BCn RSvD model fll.4p with the independent coupling 
parameters /x, ia and k. Our first goal in this section to construct a quadratic algebra relation 
for the Lax matrix A with structure coefficients satisfying the consistency conditions analogous 
to fiTTTm . 

Recalling 02.1511 and 02.8ip . we start with the observation 

BA r- 

-^^ = X2(d;a + ad:) (c6N„). 
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(4.4) 






























Therefore, upon introducing the g 0 0 -valued function 


ri 2 = 



d{h-^) 

d\c 


we can write the tensorial Poisson bracket 


{^1, ^2 — ri2-4.i -|- ^iri2. 

Now, by simply applying the Leibniz rule, from 03.104^ we get that 

{.A f* = a.i2wAi.A2 -l- *Aihi2wA2 — ^ 2 ^ 12^1 ~ *Ai.A2d!i2 


with the dynamical coefficients 

d'12 = ^^2 ^<^12^1^2 + ^^ 12 ^ 1^2 ~ ^2 ^^21^1/12, 

hl2 = ^1^2 ^^12^1 ^^2 + ^ 1 ^ 12^1 ^^2 ~ ^1^2 ^^21^.2, 
C 12 = hi ^/i 2 Ci 2 hi/i 2 ^ ^h 2 ri 2 hi -|- h2T2ihih2 
d!i2 = hih2di2hi ^/i2 ^ — /ii/i2ri2h]^ ^ -|- hih2T2ih2 ^■ 


(4.5) 


(4.6) 


(4.7) 


(4.8) 

(4.9) 

(4.10) 

(4.11) 


Since the decorations coming from h are ‘equally distributed’ among these new functions, we 
expect that likewise they satisfy the consistency conditions analogous to fid.liop . Somewhat 
surprisingly, this naive idea is fully conhrmed by the following result. 


Theorem 11. The functions ^.8 4-li ) appearing in the tensorial Poisson bracket obey 
the consistency conditions 


0'21 — —O'U, d 2 l — —di 2 , ^21 — C 12 , <2i2 + ^12 — Cl2 + <^12- 


(4.12) 


In other words, the Lax matrix A ( [.^.5] ) of the rational BCn RSvD system satisfies a quadratic 
Poisson bracket |(^. ?! ) characterized by the consistent dynamical structure coefficients (^T5[j^T7I]). 


Proof. A moment of reflection reveals that 621 = ~hi 2 , <i 2 i = ~di 2 , and ^21 = c. 12 , whence 
it is enough to prove that 0,12 -|- bu = Cu + di 2 . Since the verihcation of this last equation is 
basically an involved algebraic computation, in the following we wish to highlight only the key 
steps. First, we introduce the functions 


V{x) = \/1+'^ and Q(x) = 


{x e (0, cxd)). 

X- X 

Remembering fl4.ip . we see that 

V{x) = a{xy — I3{x)‘^ and Q{x) = 2oL{x)f3{x). 


(4.13) 


(4.14) 


To make the presentation a slightly simpler, we also introduce the G-valued function 


n = h\ 


(4.15) 
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Now, recalling fl4.ip . fl4.2p and fl4.14p . with the notations 


Va = 'P(Aa) and Qa = Q{Xa) (a e 


(4.16) 


we can clearly write 


H = 


diag(Pi, ...,Vn) diag(Qi ,..., Qr, 
-diag(Qi, ...,Qn) diag(Pi ,... ,Vn 


To proceed further, we also dehne the 0 0 0 -valued function 

^12 = -{h2'^l2 + Tuh2^)H2. 

Remembering the form of ri 2 fl4.5p . Leibniz rule yields 


C=1 

where for the derivatives we can easily find that 


n 


_idH 
dXr 


W 


:X 


26 c 


(c e Nn). 


dXc X,^yxi + K^' 

Bearing in mind the above objects, from 04.8114.11]) one can derive that 

^^2 ^(^12 + bi2 ~ Ci 2 — di2)hih2 


= Tfh ^2” “12^1^2 + ^2” ^12^2 - C12H1 - dl2 

- n^^n^2V-l + ^2~'^21^2 - f^l2 + f^21. 


(4.17) 


(4.18) 


(4.19) 


(4.20) 


(4.21) 


To handle the right hand side of the above equation, we need the commutation relations listed 
below. First, for each c G we have 


n-^D+n = T>+, (4.22) 

H-^D;n = (P| - Ql)D- + 2\PMcXi‘, (4.23) 

H-'XilH = (Vl - Ql)Xil + 2iP,Q,D:, (4.24) 

= X£l. (4.25) 

Also, if a, 6 G and a < b, then we can write 

H-'xt:i,,-H = [vj>, + s.s,)Ay;:,. + \[v,q, - nsjA'-y., ( 4 . 26 ) 

w-‘A--:.,w = (VaVi - QaQt)x:;'_,^ + i(p„Si + nsjA;y„ ( 4 . 27 ) 

n-'xi’i,^n = (v,Vt - SaSi,)A,ti„ - kp.Sj + VtOMeXc,, ( 4 . 28 ) 

w = (p.n + s.si,)A-i„ - i(p„Si - PiSj.Yiy., (4.29) 













= (p„n + QaQt,)xti„ - i(p.Sj - nejA'-i,., 
h-'x;:\-h = (p„n - s<.Si,)A-'i„ - kp.Sj + 


together with the relations 

(4.30) 

(4.31) 

(4.32) 

= (p^n + SaSi)A'-i„ + KP.Si - nSa)A«„. (4.33) 

Now, let us examine the first four terms appearing on the right hand side of fl4.2ip . Recalling 
fl3.105113.1081) . the application of 04.221) itself yields the formula 

+ 'H^' 612^2 - - di 2 

+ x-’^) A {H-^x+’^H - x+’^) 


H-^Xt:X,^H = {VaV, - QaQ,)X::l,^ + \{VaQ, + nQa)X-’:,^, 

/-I W-.i 


E 


(4.34) 


q((A) 


However, in order to further simplify this expression, we still have to exploit some functional 
equations obeyed by V and Q. By inspecting the definitions 04.13p . we see immediately that 

2 k^ 


v{xY - Q{xY = 1 + 


x^ 


(x e (0, 00 )). 


A slightly longer calculation also reveals that 

p(xYP(yY-(Q(^)Q(y)-^Y , Pi^YQiyY-nyYQi^Y 

X — y X + y 

PjxYQiyY-'PiyYQi^Y ^ H^Yny Y - iQ{x)Qiy) + 1)' 


= 0 , 


= 0 , 


(4.35) 

(4.36) 

(4.37) 


X — y X + y 

where x,y ^ (0, cxo) and x Y y- Having equipped with the relations 04.35114.37)) . let us note 
that the application of the commutation relations 04.22114.3^ does give rise to an even greater 
simplihcation in 04.34p . Indeed, we hnd that 

(«-‘A-’-W + X-'') A {H-^X+-‘H - XX‘) 


E 


q;(A) 


E 

C=1 


2 k^ 


^ 2 ;:: A Xif + 


2k 


(4.38) 


K 


^ + —Dr A X 2 Y 


A 2 


Now, let us turn to the last four terms appearing on the right hand side of 04.2ip . Recalling 
04.191) . O4.20p . and 04.351) . we can write that 


777^1212771 + 1^12 = 


C=1 


A3 


' ^ 27 : + 


2k 

a! 




' ^ 2 - 7 : 


Now, by plugging 04.34p . 04.38p and 04.39p into 04.2ip . we obtain at once that 

^^2 ^(^12 + hi2 ~ C 12 — di2)hih2 = 0, 

whence the proof is complete. 


(4.39) 

(4.40) 

□ 
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Having completed the proof, now we offer a few remarks on the resnlt. First, since the Lax 
matrix A obeys the quadratic bracket fl4.7p with the dynamical objects 04.8114.11]) satisfying the 
consistency conditions 04.12p . the quadratic bracket (\A.1\ can be rewritten as 

{^ 1 ,^ 2 }^ = [^12,*^!] ~ (4-41) 


Indeed, recalling our discussion at the end of the previous section, an appropriate r-matrix is 
provided by the formula 


ri2 = Pt2‘^2 + A 2 P 12 , 


(4.42) 


where 


, ai 2 + U12 j __ di 2 — bi 2 — C12 — U12 

Pi 2 = -^- and P 12 =- 2 - 


(4.43) 


with an arbitrary 0 V g-valued dynamical object Uu- 

Second, one may raise the objection that the formulae 04.8114.lip for the quadratic structure 
matrices in the BCn case are ‘less explicit’ than the analogous objects 03.105113.108p in the Cn 
case. The trouble is mainly caused by the derivatives of h~^ appearing in the dehnition of 
ri 2 04.5p . Though these derivatives can be worked out rather easily, we propose an alternative 
approach to cure the problem. Namely, let us apply the gauge transformation 

A{\, 9) = h{\)A{X, 9)h{X)-^ ((A, 9) e V^). (4.44) 


on the Lax matrix A 04.3p . By applying the corresponding transformation on r 04.42p . it turns 
out that the transformed r-matrix takes the form 


fi2 = P^2^2 + APi2 , 


(4.45) 


where 


= ^i^ 2 Pi 2 ^i ^^2 ^ + ^ 1^21 and = ^i^ 2 Pi 2 ^i ^^2 ^ + ^i’^ 2 r 2 i'H 2 ^ (4.46) 

Since f 12 is linear in A, the tensorial Poisson bracket for A can be cast into a quadratic form 
with structure matrices obeying the consistency conditions analogous to 04.12p . To save time 
on the algebraic details, we present only the resulting formulae. 

Theorem 12. For the transformed Lax matrix A = A'H~^ of the rational BCn RSvD model 
we have 

{.4-1, . 4 . 2 }^ = hi2.4i.42 + .4 i6i2.42 — .42Ci2.4i — .4-1.42^12, (4.47) 

where the dynamical objects 


ai2 — < 2 i 2 , 


dXr 


bi2 — Fii |^bi2 + ^77 

/ ^ n 

C 12 = 772 ( C 12 + -^ y^-p^ 


hi2=77i772 (hi2 + ^^PAA (77 


Df ) 


77 


_ih77 

dXr 


772-', 


.i<977 

dXr 


— Inj — l 




(4.48) 

(4.49) 

(4.50) 

(4.51) 
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are built up from the explicitly given functions H3.105\fJ.108\) . U-l'TD and Ii4.20\ ). Furthermore, 
by construction, they satisfy the consistency conditions 


&21 — d2l — —di2, ^21 — Ci2, &12 + &12 — ^12 + <^12- 


(4.52) 


As an immediate consequence of the r-matrix formalism, now we can easily construct a Lax 
pair for the rational BCn RSvD system. For this reason, let us recall the partial trace operation 
on the second factor of qI{N, C) 0 0 t(iV, C), which is uniquely determined by the condition 

tr 2 (X 0 F) = tr(F)X (X, F G 0l(iV, C)). (4.53) 

As we proved in |2S], for the Hamiltonian of the rational BCn RSvD model we can write 


= Itr(^) = ^tr(^). 


(4.54) 


Therefore, by expanding the Lax matrix A fl4.44|) in an arbitrary basis of the complex linear 
space 0l(iV, C), say in the basis {vj} fl2.27p . we can write 




tr(nj). 


(4.55) 


Making use of this expansion, for the action of the Hamiltonian vector held X{V^) on 

the Lax operator A we obtain 


XhAA] = J2XhAA% = 

I I 

^ I,J \ I,J 

= ^tr2 f[ri2, A] - [hi, A]) = ^[tr2(ri2), M]. 


Vi 0 vj 


(4.56) 


Since for pf 2 04.461) we can write 

^ , 0,12 + 'hl 2 1 — di2 — bi2 — Ci2 — ' 0,12 

PI 2 = -^- and P 12 =-^-, 

where 012 is a 0 V 0 -valued function, the application of 04.45p and 04.52p yields that 

tr2(f'i2) = tr2((di2 - Ci2)A). 

Upon introducing the alternative Lax matrix 

A = h-^Ah = n-^A, 


(4.57) 


(4.58) 


(4.59) 


the combination of the formulae 04.48p . 04.50^ . 04.56p and 04.58p . together with the explicit 
expressions 03.105^ and O3.107p . leads to the following result. 
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Theorem 13. With the aid of the Q-valued function 


B 


1 ^ - A))X--‘ - + A))Xi-' 

2 ^ a(X) 

Q,e ^ ^ 


1 

+ 2 tr(S'c(.4. + .4.) + Tc{A. — A.))Df 

C=1 


n 

2 


E 


tv{X£‘A) 

Ac-\/A|~+~^ 


(4.60) 


t/ie derivative of the Lax matrix A along the Hamiltonian vector field X^r takes the Lax 

form Xhr[A\ = [i3,.A]. In other words, B provides a Lax pair for A. 


5 Discussion 

One of the most important objects in the algebraic formulation of the theory of classical inte- 
grable systems is undoubtedly the r-matrix structure encoding the tensorial Poisson bracket of 
the Lax matrix. In the context of the A„-type CMS and RSvD models the underlying dynamical 
r-matrix structure is under complete control, even in the elliptic case (see e.g. [aSlEZlEHlEg]). 
In sharp contrast, for the models associated with the non-y4„-type root systems the theory is 
far less developed. By generalizing the ideas of Avan, Babelon and Talon [SI], in our earlier 
paper |22] we constructed a dynamical r-matrix structure for the most general hyperbolic BCn 
Sutherland system with three independent coupling constants. However, for the elliptic case 
only partial results are available For the non-A^-type RSvD systems the situation is even 
more delicate. Prior to our present paper, the r-matrix structure of the BCn RSvD systems 
was studied only in [H], based on the special one-parameter family of Lax matrices coming 
from Z 2 -folding of the A 2 n root system. Nevertheless, in the present paper we succeeded in con¬ 
structing a quadratic r-matrix structure for the rational BCn RSvD systems with the maximal 
number of three coupling parameters, as formulated in Theorems [TT] and [T2j It is also clear 
that by applying a standard analytic continuation argument on our formulae, one can easily 
derive a dynamical r-matrix structure for the rational RSvD system appearing in |42] . 

Regarding the hyperbolic, trigonometric and elliptic variants of the non-A„-type RSvD 
systems we also face many interesting questions. Indeed, except from some very special cases 
[2111131111], even the construction of Lax matrices for these models is a wide open problem. 
However, let us note that in the last couple of years many results for the An-tjpe models have 
been reinterpreted in a more geometric context using advanced techniques from the theory of 
reductions (see e.g. [HI 06]). Relatedly, it would be of considerable interest to see whether the 
underlying classical r-matrix structures can be explored from these geometric pictures along the 
line of our present paper. We also expect that the various reduction approaches eventually may 
lead to a progress in the rigorous geometric theory of the non-A„-type trigonometric, hyperbolic 
and elliptic RSvD systems as well. As a starting point, it is worth mentioning the recent paper 
[47] . in which a Hamiltonian reduction approach based on the Heisenberg double of SU{n,n) 
gives rise to a new integrable particle system, that in the cotangent bundle limit gives back the 
familiar hyperbolic BCn Sutherland model with three independent coupling parameters. 
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Turning back to our quadratic r-matrix algebra fl4.47|) . let us observe that the structure 
matrices 0 , 12 , 612 , C 12 and di 2 are fully dynamical, i.e. they depend on all variables of the 
phase space fll. 2 p in an essential way. It is in contrast with the CMS models, where 
the naturally appearing dynamical r-matrices usually depend only on the conhguration space 
variables. Moreover, in many variants of the CMS models the r-matrices can be related to 
the dynamical Yang-Baxter equation, as hrst realized in [48]. However, in the case Suris 
[4U] observed that in some special choice of gauge the CMS and the RSvD models can be 
characterized by the same dynamical r-matrices. Working in this gauge, Nagy, Avan and 
Rollet proved that the quadratic structure matrices of the hyperbolic An RSvD system do obey 
certain dynamical quadratic Yang-Baxter equations (see Proposition 1 in [50], and relatedly 
also [51]). As a natural next step, we hnd it an important question whether such claims can be 
made about the quadratic algebra relation fl4.47|) in an appropriate gauge. Also, it would be of 
considerable interest to investigate whether the non-A„-type RSvD models can be characterized 
by numerical, i.e. non-dynamical r-matrices. In the An case the answer is in the affirmative 
(see m, but in the BCn case the analogous tasks seem to be quite challenging even for the 
rational models. Nevertheless, we wish to come back to these problems in later publications. 
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